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ABSTRACT 

We methodically model the broad-band Suzaku spectra of a small sample of six 'bare' 
Seyfert galaxies: Ark 120, Fairall 9, MCG-02-14-009, Mrk 335, NGC 7469 and SWIFT 
J2127. 4+5654. The analysis of bare Seyferts allows a consistent and physical mod- 
elling of AGN due to a weak amount of any intrinsic warm absorption, removing the 
degeneracy between the spectral curvature due to warm absorption and the red-wing 
of the Fe K region. Through effective modelling of the broad-band spectrum and in- 
vestigating the presence of narrow neutral or ionized emission lines and reflection from 
distant material, we obtain an accurate and detailed description of the Fe K line region 
using models such as LAOR, kerrdisk and KERRCONV. 

Results suggest that ionized emission lines at 6.7 keV and 6.97 keV (particularly 
FeXXVI) are relatively common and the inclusion of these lines can greatly affect 
the parameters obtained with relativistic models i.e. spin, emissivity, inner radius of 
emission and inclination. Moderately broad components are found in all objects, but 
typically the emission originates from tens of r g , rather than within < 6 r g of the 
black hole. Results obtained with kerrdisk line profiles suggest an average emissivity 
of q^2.3 at intermediate spin values with all objects ruling out the presence of a 
maximally spinning black hole at the 90% confidence level. We also present new spin 
constraints for Mrk 335 and NGC 7469 with intermediate values of a = 0.70l^oi an d 
a = 0.69±o °g respectively. 
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1 INTRODUCTION 

Analysis of the Fe K line profile can help us determine some 
of the intrinsic properties of the central black hole in AGN 
and its accretion disc (e.g. Fabian et al. 1989; Laor 1991). 
Foremost in the current climate is the determination of black 
hole (BH) spin, for example recent observations by Miniutti 
et al. (2007 & 2009) and Schmoll et al. (2009) have analysed 
the broad Fe Ka region in order to constrain the spin of the 
central BH. A black hole can be characterised simply by its 
mass and its spin. Many objects have now been classified in 
terms of their mass into three categories: Galactic BH, In- 
termediate Mass BH and Supermassive BH. The spin of the 
black hole determines the nature of the space-time metric 
in the regions close to it. The spin parameter a = cJ/GM 2 
(where J=angular momentum and < a < 0.998) is used 
to describe the spin of the BH, where a is limited to a max- 
imum value of a = 0.998 at the Thorne limit. This is due 
to photon capture in which photons travelling on 'negative' 
angular momentum orbits are preferentially captured by the 



black hole therefore producing an upper bound to the spin 
parameter a and hence limiting the innermost stable circu- 
lar orbit to a minimum of rrsco = 1.235 r g (Thorne 1974). 
Constraining the spin of SMBHs in AGN and studying the 
distribution of black hole spin can aid our understanding of 
the evolution of AGN and the black holes themselves e.g. 
mergers, relativistic jets and variability (Blandford & Zna- 
jek 1977; Volonteri et al. 2007; King, Pringle & Hofmann 
2008). 

Line emission from the inner regions of the accretion 
disc can become broadened due to relativistic effects and 
Doppler motions (Fabian et al. 1989) resulting in an asym- 
metric profile. Evidence for such broadening was found using 
the X-ray CCD detectors onboard ASCA (Tanaka, Inoue & 
Holt 1994) by Mushotzky et al. (1995), Tanaka et al. (1995) 
and Nandra et al. (1997) typically over the 0.5-10.0 keV 
range. Now with X-ray spectra of increasing quality and over 
wider energy ranges such as that obtained with the XIS and 
high energy HXD detectors onboard Suzaku (0.5-60.0 keV, 
Mitsuda et al. 2007), the FeK region of AGN can be exam- 
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Table 1. The Suzaku Scyfcrt sample 



Object 


RA (J2000) 


Dec (J2000) 


Redsbift 


N H (Gal) (10 22 cm- 2 ) 


Ark 120 


05 16 11.4 


-00 09 00 


0.033 


0.0978 


Fairall 9 


01 23 45.8 


-58 48 21 


0.047 


0.0316 


MCG-02-14-009 


05 16 21.2 


-10 33 41 


0.028 


0.0924 


Mrk 335 


00 06 19.5 


+20 12 10 


0.026 


0.0356 


NGC 7469 


23 00 44.4 


+08 36 17 


0.016 


0.0445 


SWIFT J2127.4+5654 


21 27 45.0 


+56 56 40 


0.014 


0.7650 



Table 2. Summary of observations for the objects in the sample. 1 The observed 2-10 keV flux for XIS and EPIC-pn instruments, 
15-50 keV flux for HXD and 20-100 keV flux for BAT, in units lO^ci-gcm" 2 s" 1 from Model A. 



Object 


Mission 


Instrument 


Date 


Exposure (s) 


Count rate 




Flux 1 


Obs. ID 


Ark 120 


SuzaJcu 


XIS 
HXD 




2007/04/01 


100864 
89470 


1.896 ± 0.003 
0.114 ± 0.003 




3.06 
3.46 


702014010 


XMM 


EPIC- 


-pn 


2003/08/24 


78170 


25.04 + 0.018 




3.87 


014719010 




Swift 


BAT 






2453000 


(6.9 + 0.4) x 10~ 


-4 


4.88 




Fairall 9 

i_ CAtll CXI 1 c/ 


SuzaJcu 


XIS 
HXD 




2007/06/07 


167814 
127310 


1.718 ± 0.002 
0.089 ± 0.002 




2.53 
2.97 


702043010 


XMM 


EPIC- 


-pn 


2000/07/05 


25830 


5.825 + 0.015 




1.17 


0101040201 




Swift 


BAT 






3280000 


(4.9 + 0.4) x 10~ 


-4 


4.35 




MCG-02-14-009 


SuzaJcu 


XIS 
HXD 




2008/08/28 


142152 
120028 


0.216 + 0.001 
0.017 + 0.002 




0.42 
0.61 


703060010 




XMM 


EPIC- 


-pn 


2009/02/27 


82790 


1.946 + 0.005 




0.43 


0550640101 


Mrk 335 


SuzaJcu 


XIS 
HXD 




2006/08/21 


151296 
131744 


1.324 + 0.002 
0.012 + 0.001 




1.49 
1.31 


701031010 


XMM 


EPIC- 


-pn 


2006/01/03 


80360 


16.43 + 0.014 




1.82 


0306870101 




Swift 


BAT 






3273000 


(2.5 + 0.3) x 10~ 


-4 


1.74 






SuzaJcu 


XIS 
HXD 




2008/06/24 


112113 

85315 


1.091 + 0.002 
0.068 + 0.002 




2.11 
3.22 


703028010 


NGC 7469 


XMM 


EPIC- 


-pn 


2004/11/30 


53014 


24.09 + 0.021 




2.74 


0207090101 




EPIC- 


-pn 


2004/12/03 


54977 


16.32 + 0.017 




2.89 


0207090201 




Swift 


BAT 






3286000 


(6.6 + 0.3) x 10- 


-4 


4.82 




SWIFT J2127.4+5654 


SuzaJcu 


XIS 
HXD 




2007/12/09 


91730 
83321 


1.373 + 0.004 
0.074 ± 0.002 




3.77 
3.33 


702122010 




Swift 


BAT 






3999000 


(4.2 + 0.3) x 10- 


-4 


4.16 





ined in detail. High energy X-ray data are important since it 
allows the reflection component and its strength to be prop- 
erly fit, assessing its contribution to the continuum and Fe K 
region (e.g. Reeves et al. 2007). Fitting features such as the 
Compton hump at ~ 30keV allows, for example, the ioniza- 
tion state of the reflecting material to be determined (Ross 
& Fabian 2005). With the aim of measuring properties of the 
accretion disc and the central black hole itself, broad-band 
data allows us to start making constraints on parameters in 
these regions based upon the shape of the Fe K line profile. 

Previous studies of iron lines have been made using 
data from XMM-Newton over the 2.5-10.0 keV (Nandra et 
al. 2007) and 0.6-10.0 keV energy ranges (Brenneman & 
Reynolds 2009) , finding complex emission in the Fe K band 
in the majority of Type 1 Seyfert AGN over and above nar- 
row line components originating from distant material. In a 
sample of 26 objects Nandra et al. (2007) found that nar- 
row 6.4 keV emission is ubiquitous amongst AGN and broad 
FeK lines feature in ~40-50% of AGN and ionized emis- 
sion due to Fe XXVI and Fe XXV is relatively rare amongst 



AGN. Brenneman & Reynolds (2009) found that 4/8 AGN 
were best fit by a model consisting of relativistically blurred 
reflection from the inner regions of the accretion disc with 
2/8 objects suggesting non-zero spin, however noting that 
the 10.0 keV cutoff with the EPIC-pn camera limits their 
findings of the reflection continuum. 

This small sample features predominantly bare Seyferts 
i.e. with a very weak warm absorber. This is in an attempt 
to simplify the modelling of the broad-band spectrum and 
therefore provide a basic understanding of the properties 
of the FeK region and any accompanying broad red-wing 
without the need for debate over differing interpretations of 
the origin of various absorption components in the spec- 
trum (Turner & Miller 2009). Using data from Suzaku's 
XIS (Koyama et al. 2007) and HXD (Takahashi et al. 2007) 
detectors spanning 0.5-60.0 keV and BAT data (from the 
Swift 22 month all sky survey, Tueller et al. 2010) over 
20.0-100.0 keV provides the broad-band spectra necessary 
for detailed modelling and measurement of the Fe K region 
and the associated Compton reflection hump. This approach 
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ensures that a robust model of the FeK line region (and 
broad-band spectrum) can be made with the aim of estab- 
lishing the degree to which narrow ionized emission lines and 
relativistically broadened components are required. 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Observations 

The objects in this sample are described in Table [T] and are 
all the Seyfert 1, radio quiet AGN with low intrinsic ab- 
sorption available in the public Suzaku data archive with 
spectra from the XIS and HXD detectors. All objects are 
detected above 15 keV and are relatively nearby with red- 
shift 2 < 0.05. Hard X-ray data from the BAT instrument 
onboard Swift is used in addition to the HXD data in all 
observations other than MCG-02-14-009, for which no BAT 
spectrum was publically available (it was however detected 
by BAT in the 39 month survey with a 14-150 keV flux of 
1.5 x 10 -11 erg cm -2 s _1 , Cusumano et al. 2010). In addition 
to the Suzaku data, observations with XMM-Newton are 
also used for comparison, except for SWIFT J2127.4+5654 
which was not publically available. The details of the obser- 
vations used in this analysis are given in Tabled 

The objects featuring in this sample, however, are all 
objects with a very low degree of warm absoprtion. Three 
objects in this analysis have been noted as having some addi- 
tional absorption in previous observations: MCG-02-14-009, 
Mrk 335 and NGC 7469. Gallo et al. (2006) found a small 
neutral absorber (Ah < 10 21 cm -2 ) and a OVII edge in 
a re-analysis of a short (~ 5 ks net exposure) observation 
of MCG-02-14-009 with XMM-Newton. However including 
these components in the Suzaku data used here makes no 
improvement (Ah < 3 x 10 20 cm -2 ) and the optical depth of 
the OVII edge can only be constrained to r < 0.04. There- 
fore MCG-02-14-009 can be considered as 'bare' for the pur- 
poses of this paper. 

Mrk 335 has also shown evidence for a warm emitter 
whilst in a low state when observed with XMM-Newton 
(Grupe et al. 2008), however in the Suzaku data used here 
Mrk 335 has a 0.5 - 2.0 keV flux 14x higher than in the low 
state observation and as such any warm emission features 
are entirely dominated by the continuum. As a result of this 
Mrk 335 is also suitable for inclusion within this sample. 
Previous observations of NGC 7469 have noted some degree 
of X-ray and UV absorption (A H ~ 10 20 - 10 21 cm -2 , Scott 
et al. 2005; Blustin et al. 2007). The Suzaku data used here 
is consistent with the previous work using data from XMM- 
Newton, requiring an OVII edge depth of r < 0.1, however 
there is no effect upon the FeK parameters and NGC 7469 
has also been included within this sample. Figure[5]shows no 
significant additional absorption features below 2 keV within 
the Suzaku data. 



2.2 Data reduction 

All the Suzaku data in this paper were reduced using the 
HEASOFT reduction and analysis package (version 6.8). 
XIS source spectra were extracted from circular regions of 
3.0' within XSELECT centred upon the source at the on- 
axis pointing position. Similarly background spectra were 



extracted from 3.0'circular regions, taking care not to in- 
clude the source or the Fe55 calibration sources in the 
corners of the CCD's field of view. Only data from the 
front-illuminated XIS cameras were used i.e. the XIS0 and 
XIS 3 cameras - due to their greater sensitivity at FeK 
energies, however the observation of Mrk 335 (Obs. ID 
701031010) also includes the now non-operational front- 
illuminated XIS 2. It should be noted that only data from 
the XIS 3 camera was available for the Suzaku observation of 
SWIFT J2127.4+5654 (Obs. ID 702122010) since the XIS 
camera was not operating properly. 

XIS redistribution matrix files (rmf) were created using 
the HEASOFT tool xisrmfgen and the ancilliary response 
files (arf) using xissimarfgen. Data from the 3x3 and 5x5 
modes were grouped together and the data from each of the 
XIS front-illuminated CCDs were co-added in order to in- 
crease signal to noise using mathpha, addrmf and addarf. 
We ignore all XIS data below 0.5 keV, above 10.0 keV and 
between 1.7-1.95 keV due to uncertainties in the calibration 
of the detectors around the SiK edge. 

The HXD/PIN spectrum was extracted from the 
cleaned HXD/PIN events files and corrected for dead time 
using the tool hxddtcor. The tuned HXD/PIN background 
events were used for background subtraction (Fukazawa et 
al. 2009) using identical good time intervals (GTIs) as per 
the source events, generated with 10 x the actual back- 
ground count rate to minimize photon noise. A simulated 
cosmic X-ray background was also produced using XSPEC v 
12.5. In with a spectral form identical to Gruber et al. (1999) 
and then added to the corrected non X-ray background file 
to create a single background file. In the analysis of the 
HXD/PIN spectra we typically consider data between 15.0- 
60.0 keV however in some cases such as Mrk 335 the data 
was also ignored below 20.0 keV due to thermal noise. Hard 
X-ray data also obtained from the Swift 22-month BAT cat- 
alogue was included for all objects (other than MCG-02-14- 
009) over the energy range 20.0-100.0 keV. 



3 ANALYSIS & RESULTS 

Spectral analysis and model fitting is performed with 
XSPEC v 12.5. In (Arnaud 1996), all models were modified 
by Galactic absorption via the WABS multiplicative model 
(Morrison & McCammon 1983) using a Galactic column 
density obtained using the nh ftool appropriate for each 
source giving the weighted average Ah value of the (LAB) 
Survey of Galactic HI (Kalberla et al. 2005), using abun- 
dances from Anders & Grevesse (1989). In all fits a constant 
factor was introduced to account for the cross-normalization 
between the XIS, HXD/PIN and BAT detectors; fixed at 
1.16 or 1.18 between XIS and HXD/PIN according to the 
nominal pointing positiorQ and allowed to vary between XIS 
and BAT detectors (typically ~ 1.0, indicating little vari- 
ability between the Suzaku and BAT datasets). Data is fit 
over the full 0.5-100.0 keV range, excluding those regions 
mentioned above. The % 2 minimization technique is used 
throughout, quoting 90% errors for one interesting parame- 
ter (Ax 2 = 2.71) unless otherwise stated. Where the signifi- 

1 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/watchout.html 
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Figure 1. The 0.5-70.0 keV spectrum and residuals after modelling of the continuum with a simple POWERLAW and WABS to account for 
Galactic absorption, the entire FeK region and any soft excess is left unmodelled. XIS data is in black, HXD in red and BAT data is 
represented by blue circles. 



Table 3. Summary of models and components used for each feature. 



Model 


Continuum 


Distant Reflection 


Soft Excess 


6.4 keV Core 


6.4 keV Broad 


Model A 


POWERLAW 


PEXRAV 


COMPTT 


Gaussian 


Gaussian 


Model B 


POWERLAW 


REFLIONX 


compTT 


within REFLIONX 


None 


Model C 


POWERLAW 


REFLIONX 


compTT 


within REFLIONX 


LAOR 


Model D 


POWERLAW 


REFLIONX 


compTT 


within REFLIONX 


KERRDISK 


Model E 


POWERLAW 


REFLIONX 


Blurred REFLIONX 


within REFLIONX 


Blurred REFLIONX 


Model F 


POWERLAW 




compTT 


Gaussian 


Blurred REFLIONX 



cance of components is stated according to A\ 2 , the compo- 
nent has been subtracted from the final model and then refit 
to ensure the order in which components are added does not 
effect the quoted statistcal significance. Table [3] summarises 
the subsequent models and the components used to model 
the continuum, soft excess and FeK region. 



3.1 Model A — Parametrization of spectra 

Model A is intended to provide a simple parametrization 
of the spectra and give an insight into the presence of ba- 
sic components of the spectra and the extent to which a 
possibly relativistically broadened component is required to 
model the Fe K line region. None of the spectra required any 



© 2010 RAS, MNRAS 000, ??-?? 



Iron line profiles in Suzaku spectra of bare Seyfert galaxies 5 




1.4 
1.2 
1 

0.8 









1 10 

Observed Energy (keV) 
1 1 1 i i, i 1 iii 



10 

Observed Energy (keV) 




1.4 
1.2 
1 

0.8 



Observed Energy (keV) 



Observed Energy (keV) 



Figure 2. The 0.5-70.0 keV spectrum and residuals after modelling of the continuum with a POWERLAW, neutral reflection from the 
PEXRAV model, COMPTT to model the soft excess and WABS to account for Galactic absorption, the entire Fe K region is left unmodelled. 
XIS data is in black, HXD in red and BAT data is represented by blue circles. 



significant warm absorber (see Section 2.1), thereby simpli- 
fying any present broadening in the region. Additionally, the 
COMPTT model representing Comptonization of soft pho- 
tons in a hot plasma above the disc (Titarchuk 1994) with 
a soft photon input temperature of 0.02 keV, is employed to 
account for the soft excess if present in the spectra (see Fig- 
ure [T}. Porquet et al. (2004) also found that in a sample of 
PG quasars the soft excess was better modelled in this way, 
rather than by thermal emission from the accretion disc. 
A second soft POWERLAW component instead of COMPTT 
also gives a similar parametrization of the soft excess. The 
narrow 6.4 keV core due to reflection from distant material 
is present in all six objects and has been modelled with a 
narrow Gaussian with width r/Ko free to vary. The narrow 
component is not resolved in any of the spectra and as such 
the width is fixed at r/Ko = 0.01 keV. Emission resulting 
from FeK/3 is also accounted for with the line energy fixed 



at 7.056 keV, width fixed to that of the narrow Ka and flux 
tied to 13% of the Ka component. 

Consistent with this, neutral distant reflection is ac- 
counted for using the pexrav model (Magdziarz & Zdziarski 
1995) applied to the broad-band 0.5-100 keV spectra. This 
model requires the input of a photon index T which is 
tied to the continuum powerlaw, the normalization of the 
pexrav component is also tied to that of the powerlaw, 
abundances are assumed to be Solar (Anders & Grevesse 
1989) and the disc inclination to the observer is fixed at 
cosi = 0.87 throughout. The reflection fraction R = S1/27T 
is left as a free parameter (where R = 1 denotes reflection 
from material subtending 2-7rsr). The cut-off energy for the 
pexrav component is fixed at 1000 keV, except for NGC 
7469 and SWIFT J2127.4+5654 which show indications of a 
roll-over at high energies, occuring at E c — 1191^ keV and 
E c = 491*4 keV respectively. Such a roll-over was found in 
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Figure 3. The 4-9 keV residuals after modelling of the continuum with a powerlaw, COMPTT to model the soft excess and Galactic 
photoelectric absorption. The dashed vertical line represents 6.4 keV in the rest frame. 



SWIFT J2127.4+5654 previously by Malizia et al. (2008) at 
E c = 33±H keV. 

Examining the residuals after the application of these 
components, some objects have more complex features 
such as excess emission around energies of 6.7 keV and 
6.97 keV relating to narrow ionized emission from FeXXV 
and Fe XXVI respectively, again due to distant photoionised 
gas (see Figure [2] and Figure [3} . Similarly to the modelling 
of the narrow 6.4 keV core, the energies of these lines are free 
to vary as is the normalization, however the width remains 
fixed to a = 0.01 keV. A broad Gaussian is also added to the 
model to account for a red-wing in the spectra with energy, 
normalization and rJBroad as free parameters. In this case the 
energy is allowed to drop below 6.4 keV in the rest frame as 
this model is only intended as a simple parametrization of 
the Fe K region and a test of the significance of the compo- 
nents used to model the region. 

A good fit is obtained by Model A in all objects, par- 
ticularly Fairall 9, NGC 7469 and SWIFT J2127.4+5654. 



No significant soft excess is found in MCG-02- 14-009 and 
SWIFT J2127.4+5654, the latter (if present) is likely to 
be mostly absorbed due to the relatively high amount of 
Galactic absorption (jV H ~ 7.65 x 10 21 cm _2 ). In SWIFT 
J2127. 4+5654 there is some indication of the presence of an 
intrinsic neutral absorber, albeit of relatively low column 
density (as found in Miniutti et al. 2009). Instead mod- 
elling the spectrum with an absorbed powerlaw using the 
zphabs model at the redshift of the source improves the fit 
by Ax 2 ~ 59 for one addition free parameter with intrinsic 
column density Nh ~ 8.2l{j'g x 10 20 cm -2 and photon index 

Statistically significant narrow ionized emission is found 
in most objects, with SWIFT J2127.4+5654 and MCG- 
02- 14-009 having very little requirement for these features. 
These lines occur at energies likely originating from Fe XXV 
and FeXXVI with emission at ~ 6.97 keV being particu- 
larly evident in the spectra, proving strong in Fairall 9 for 
which the fit improves by A^ 2 ~ 24 with the introduction 
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Table 4. Model A components for Suzaku XIS, HXD and BAT data from Swift. a POWERLAW normalization given in units 
(1(T 2 phkeV" 1 cm" 2 s" 1 ). b Flux for COMPTT quoted over the 0.5-10.0 keV range in units 10 _11 crgcm" 2 s _1 . 
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of a narrow Gaussian (and therefore two additional free pa- 
rameters) at 6.981q'q2 keV and to a lesser extent in Ark 120 
(Ax 2 ~ 12 at 6.96+ ;^keV). The only objects not show- 
ing residuals at ~ 6.7 keV or ~ 6.97 keV are NGC 7469 and 
SWIFT J2127.4+5654, also in agreement with Miniutti et 
al. (2009). Residuals at ~ 6.7 keV are only found to be sig- 
nificant in two objects: Fairall 9 and Mrk 335 with improve- 
ments of Ax 2 ~ 7 and ~ 14 respectively. 

With the aim of determining the extent to which 
emission from further in to the black hole is required, a 
broad component significantly improves the quality of the 
fit for most objects, typically Ax 2 > 20 for three addi- 
tional free parameters and with line widths of the order 
o-Broad >0.3keV. Only in MCG-02-14-009 does the addition 
of a broad Gaussian not make a particularly significant im- 
provement (Ax 2 ~ 3). This is surprising given evidence 
to the contrary by Porquet (2006) in which a prominent 
broad and statistically significant iron line was found in a 
short XMM-Newton observation with an equivalent width 



of EW~ 527+JIs eV whereas only EW~ 92±|| eV is found 
here. However the XMM-Newton observation was only 5ks 
net exposure and as a result the FeK line parameters are 
poorly constrained, while the Suzaku data also allow the 
broad-band continuum to be better constrained (however 
we cannot rule out some variability between these two ob- 
servations). In SWIFT J2127.4+5654, this may be due to 
the relatively high best-fitting value of the reflection compo- 
nent with R = 2.5l ' 8 which may reduce the significance of a 
broad component. Mrk 335 features a relatively broad Gaus- 
sian with tTBroad = 0.50+.°,;" keV and EW= lU±f s eV, how- 
ever this feature is not as strong as the one found by Larsson 
et al. (2008) with a = 0.45±g;oe keV and EW= 250t^ eV. 
Larsson et al. (2008) also find that the FeXXVI emission 
line does not improve the fit, however an improvement of 
Ax 2 ~ 8 for three additional free parameters is found here 
with the introduction of a line at 6.98lo'i4 keV. 

In general, the broad emission found in these objects 
is typically too strong to be modelled as purely a Comp- 
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Table 5. Model B components for Suzaku XIS, HXD and BAT data from Swift. The ionization parameter £ is given in units erg cms 1 . 
a POWERLAW normalization given in units (10 — 2 phkeV -1 cm -2 s _1 ). b REFLIONX normalisation given in units 10~ 5 . 
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ton shoulder despite the similarity of the line energy of the 
broad Gaussian with the first order Compton shoulder peak 
energy of ~ 6.24 keV (see Table [4}. However we cannot rule 
out the possibility of a contribution from a Compton shoul- 
der. Given the relatively high EW of the broad line com- 
ponent (mean EW ~ 86 eV) compared to the narrow core, 
the majority of the emission is likely to arise from a broad- 
ened component. Indeed the Compton shoulder is unlikely 
to contribute more than ~ 20% of the 6.4 keV core flux (e.g. 
George & Fabian 1991; Matt 2002). Note that the Compton 
shoulder is included in subsequent models through the use 
of a reflionx reflection component. 



3.2 Model B — Self consistent distant reflection 

In an attempt to build a more self-consistent model the 
pexrav and narrow 6.4 keV core are replaced with the re- 
flionx ionized reflection model from Ross & Fabian (2005) 
which includes emission lines in addition to the reflection 
continuum. Model B provides fits of a reasonable quality 
to all six spectra, however it is not expected to improve 
upon Model A since it doesn't include a broad Gaussian 
to parametrize the red-wing in an attempt to retain self- 
consistency and begin the construction of a more physically 
motivated model. Model B therefore acts as a 'null hypothe- 
sis' model in that all emission and reflection originates from 
distant matter. In most cases the fit statistic is somewhat 
worse than in Model A, except for Fairall 9, in which the 
quality of the fit is actually improved to xl = 864.7/834. 
The COMpTT model is still used here to model the soft ex- 
cess in the four objects requiring it, however the COMPTT 
parameters are not quoted in Table [S] (similarly in future 
tables) since the parameters obtained are consistent with 
those in Model A. 



The narrow ionized emission lines at 6.7keV and 
6.97 keV found in Model A remain in Model B and following 
the removal of the broad component the significance of these 
lines is increased. In the case of SWIFT J2127.4+5654, for 
which there are no statistically significant ionized emission 
lines in Model A, there are small excesses at these energies, 
indeed suggesting the possible presence of these lines. Emis- 
sion lines of high significance in Model A are also similarly 
significant in Model B e.g. the FeXXVI in Fairall 9 is par- 
ticularly prominent in both A & B. This is contrary to the 
analysis of Fairall 9 by Schmoll et al. (2009) in which they 
did not detect the FeXXVI line in the same Suzaku data, 
but they did however detect the He- like FeXXV emission 
line which is also seen here in Models A and B. 

All objects show some indication of at least a small red- 
wing (see Figure g}, noting that in MCG-02-14-009 this fea- 
ture is not particularly strong and in all objects the red wing 
does not extend below ~ 5 keV. A smaller than expected ex- 
cess at energies red- ward of 6.4 keV is seen in Fairall 9, given 
the statistical significance of the broad Gaussian in Model A. 
This could be due in part to a relatively high Fe abundance 
as measured by the Fe/Solar parameter within reflionx 
at approximately twice the solar value. Some objects in the 
sample indicate a super-solar Fe abundance according to the 
reflionx component (see Table O, however the true abun- 
dance may in fact be lower than this measured value if part 
of the Fe Kq emission additionally arises from Compton-thin 
matter such as the BLR or NLR. A summary of the results 
obtained with Model B can be found in Table [5] 

3.3 Model C - Laor profile 

The presence of emission from close-in to the black hole has 
been indicated in the previous two models. In accordance 
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Figure 4. Ratio plots of Model B (i.e. without including a broad FeKct) revealing excesses at energies red- ward of 6.4 keV in some 
objects. Model B consisting of POWERLAW + COMPTT (where required) + unblurred REFLIONX + narrow ionized lines as required. 



with this, emission from the FeK line region can be mod- 
elled with a relativistic component representing line emis- 
sion from the accretion disc in addition to any appropriate 
emission from distant material (see Figure |4}. Model C ad- 
vances upon Model B by remodelling the FeK region with 
the addition of a LAOR line profile operating under the as- 
sumption of a maximally spinning black hole (Laor 1991). 
The line energy was restricted to 6.4-6.97 keV in the rest 
frame, with emissivity, inclination and the inner radius of 
emission allowed to vary. The outer radius of emission was 
fixed at 400 r g throughout. In the cases where the line emis- 
sion reaches its lower limit it is fixed at 6.4 keV in the rest 
frame. Whilst giving a more feasible interpretation of the 
broad emission in the FeK region it does not provide us 
with the most physically accurate representation given the 
assumption of a maximally rotating central black hole with 
a — 0.998 within the LAOR model. Not all of the objects 
in the sample will feature a maximally rotating black hole, 
indeed Fairall 9 and SWIFT J2127.4+5654 have previously 



been found to have intermediate spin values (see Schmoll et 
al. 2009; Miniutti et al. 2009). The presence of narrow ion- 
ized emission lines was reassessed after accounting for the 
blue- wing of the LAOR profile. Employment of this model 
therefore seeks to further parametrize the broad emission in 
the Fe K region, providing suitable and plausible parameters 
for use in the later kerrdisk models, whilst giving us an in- 
dication of the extent to which the spin of the central black 
hole has an effect upon the observed spectrum. 

The fit to all objects is improved over the purely distant 
emission in Model B. The 6.97 keV line is found to be present 
in all but NGC 7469 and SWIFT J2127.4+5654 whilst the 
6.7 keV line is only found in Mrk 335 (although relatively 
weak). Model C improves the fit to MCG-02- 14-009 the least 
with only A\ 2 ~ 14 for three additional free parameters 
whereas the introduction of a LAOR profile offers a significant 
improvement for most other objects, particularly Fairall 9, 
Mrk 335 and SWIFT J2127.4+5654 (Ax 2 ~ 33, 37 & 34 
respectively), see Table [6] 
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Table 6. Fit parameters from Model C to Suzaku XIS, HXD and BAT data from Swift. Line energies are quoted in the rest frame. * 
denotes a frozen parameter, for cases where the emissivity index is fixed at q = 3 it is unconstrained. The improvement Ax 2 in the fit 
with the introduction of a LAOR profile is noted in comparison with the purely distant reflection as present in this model. Note that the 
FcXXV emission line is no longer required in Ark 120, Fairall 9 and SWIFT J2127. 4+5654. a POWERLAW normalization given in units 
(10 -2 phkcV -1 enr 2 s _1 ). b REFLIONX normalisation given in units 10~ 5 . 
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None of the objects in the sample require emission from 
within 6 GM/c 2 indicating that a rotating black hole, whilst 
possible, is not required to model the spectra. According to 
the fit parameters obtained with LAOR, the inner radius of 
emission lies at tens of r g for all six AGN. In most cases 
the accretion disc is unlikely to be truncated at these dis- 
tances from the black hole and it is likely that these values 
arise from the assumption of a maximally rotating black 
hole within the LAOR model. The emissivity indicies also in- 
dicate that a high concentration of emission from very close 
to the black hole is absent in the spectra from these AGN. 
An emissivity index of q > 5 would suggest that the emis- 
sion from the accretion disc is very centrally concentrated 
(e.g. Miniutti et al. 2003), while for objects such as Mrk 
335 and Fairall 9 this is within the error bars, one would 
also expect the inner radius of emission to be much closer 
to the black hole if this interpretation were suitable for the 
spectra. For the cases where the emissivity index cannot be 
constrained, it is fixed to q = 3.0 which is consistent with 
the other objects in the sample. 

Mrk 335 appears to have one of the strongest relativis- 
tic lines in the sample with a LAOR profile equivalent width 
of EW = I26I35 eV . and an improvement of Ax 2 ~ 37. 



However this is not particularly strong in comparison with 
previous studies of this AGN, for example Longinotti et 
al. (2007) find an EW = 320±\™ eV at a line energy of 
E = 6.9318:27 keV in a 40 ks XMM-Newton observation of 
Mrk 335 in 2000, also using a LAOR profile. The width of 
the iron line emission from the disc and its suggested high 
ionization state could be due to the lack of a Fe XXVI line in 
the Longinotti spectra (possibly due to low S/N), which has 
been accounted for in this analysis and is noted by O'Neill et 
al. (2007) in an analysis of the same XMM-Newton observa- 
tion used here. Also in agreement with our results, O'Neill et 
al. found the equivalent width of the broad line in Mrk 335 to 
be EW = 1151" eV in comparison with EW = 126±|| eV 
for the LAOR profile and EW = H3I52 eV for the broad 
Gaussian employed in Model A. 

The results obtained here are consistent with a recent 
analysis of Ark 120 by Nardini et al. (2010, submitted) who 
initially fit the features in the Fe K region with a LAOR pro- 
file. In the Nardini et al. (2010) analysis, both the inclination 
of the accretion disk and emissivity are frozen at typical 
values of i — 40° and q — 3.0 respectively, finding an in- 
ner radius of emission of r- ln = 13l 7 9 r g . These values are 
consistent with those obtained here with Model C, finding 
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i° = 34+1, q = 3.0±l l and r in = 25±J, 9 r g . Also detected is 
the FeXXVI narrow ionized emission line, again consistent 
with Nardini et al. (2010). 

As in Model A, there are no narrow ionized emis- 
sion lines found in SWIFT J2127. 4+5654 corresponding to 
FeXXV and FeXXVI. This implies that the excesses mod- 
elled as narrow components in Model B may instead be due 
to relativistically broadened FeKa emission. When mod- 
elled as a single LAOR profile all significant excesses in the 
Fe K band are removed with the profile centroid rest frame 
energy at 6.4 keV and the blue- wing peaking at ~ 6.7 keV in 
the rest frame. Consequently the observed excess at 6.7 keV 
observed in Model A and more significantly in Model B is 
more convincingly described as a blue-wing. Additionally, 
the observed intermediate feature in MCG-02- 14-009 is mod- 
elled well by the blue-wing of the LAOR profile which occurs 
at ~ 6.58 keV, prompting the removal of the previously em- 
ployed narrow Gaussian. 

3.4 Model D - Kerrdisk profile 

Replacing the LAOR model with a kerrdisk (Brenneman & 
Reynolds 2006) line profile further allows a more physically 
motivated fit to the data. The kerrdisk model allows the 
spin parameter to be varied between < a < 0.998 with the 
aim of determining the extent to which we can rule in or out 
a non-rotating or maximally spinning black hole. Similarly 
to Model C, the outer radius of the disc is fixed at 400 r ms , 
while assuming a disc of uniform emissivity. Throughout the 
fits using Model D the inner radius of emission is assumed to 
extend down to the innermost stable circular orbit (risco)- 
The line energy is confined to 6.4-6.97 keV in the rest frame, 
being frozen at 6.4 keV if it reaches its lower limit. 

Model D provides good fits to all six objects, producing 
the best fitting physically motivated model for Mrk 335 and 
NGC 7469 prior to considering a blurred reflection compo- 
nent from the inner regions of the accretion disc (i.e. Models 
E & F), see Table [7] and Figure [5] The presence of ionized 
emission lines is entirely consistent with Model C, adding 
weight to the probability that these lines are present in the 
spectra. Further to this, the inclinations of the accretion 
discs to the observer are all comparable to those obtained 
in Model C, as are the equivalent widths of the relativistic 
line profiles. Given that the spin of the central black hole 
is a free parameter and the assumptions within this model, 
the emissivity indicies are expected to vary from those ob- 
tained with the LAOR profile since the measured emissivity 
is degenerate to some extent with the inner radius of emis- 
sion. Therefore as the spin parameter varies so does risco, 
consequently affecting the measured q. 

The results obtained from the spin parameter a for these 
six objects suggest that a maximally spinning central black 
hole can be ruled out at the 90% confidence level, however 
for some objects such as Ark 120 and MCG-02-14-009 the 
spin is unconstrained and only an upper limit can be mea- 
sured (a < 0.94 and a < 0.88, suggesting that emission does 
not occur within 2.02 r s and 2.45 r s respectively, consistent 
with the Nardini et al. 2010 analysis of Ark 120) . 

Here for Fairall 9 the best fitting model gives a measured 
spin value of a — 0.441q'^. A previous spin constraint for 
Fairall 9 by Schmoll et al. (2009) found a = 0.65±2;S1 using a 
blurred reflector model (kerrconv, Brenneman & Reynolds 



(2006), convolved with reflionx). The measurement here 
also gives an intermediate value, however it is only consis- 
tent with their findings when they ignore the spectra below 
2keV to ensure that the soft excess is not the component 
driving the main part of the fit. Schmoll et al. (2009) quote 
a worse spin constraint when these conditions are upheld 
a — O.5I03. Within error bars at the 90% level these re- 
sults are consistent with the value of the spin parameter 
found here. However, here we find the emissivity index is 
constrained to q — 2.7+q'I whereas Schmoll et al. only con- 
strain this to q > 4.9. 

Mrk 335 has been noted previously as an object with a 
broad relativistic FeK line (Gondoin et al. 2002). Here we 
measure an intermediate spin a — 0.70^0 qi, ruling out max- 
imally spinning and non-rotating black holes only at the 95% 
confidence level. The emissivity index is measured at a mod- 
erate value of q = 2.6^Q3 in agreement with previous stud- 
ies of this object, for example by Longinotti et al. (2007). 
In agreement with Model C, the kerrdisk component also 
features a relatively broad equivalent width EW = 146lgg 
with an improvement of A\ 2 ~ 40 (for four additional free 
parameters) over purely distant emission. 

In NGC 7469 we measure a broad Fe line featuring a low 
emissivity index of q — 1.71q g, an inclination of i° — 23 + _\ 5 
and EW = 9li| eV, providing good agreement with the 
previous measurements made within Model C. We obtain a 
black hole spin of a = 0.691qq9, ruling out a Schwarzschild 
black hole with 95% confidence and a maximally rotating 
black hole at the 99% confidence level. 

SWIFT J2127.4+5654 also shows a broad and statis- 
tically significant component, EW = 178^gg eV, improving 
the fit by A\ 2 ~ 37 (similarly to Model C) for four addi- 
tional free parameters. The emissivity and the inclination of 
the accretion disc are typical of the objects in this sample. 
The spin parameter measured here is consistent with that 
measured by Miniutti et al. (2009), 0.70±g;™ 

compared to 

0.6ltQ 2 found using a blurred reflection model. Here we re- 
ject a maximally rotating central black hole at greater than 
99% confidence and a non-rotating black hole with 95% con- 
fidence. The spin constraints obtained with Model D for Mrk 
335, NGC 7469 and SWIFT J2127.4+5654 can be seen in 
Figure [7] 



3.5 Model E - Blurred reflector (1) 

In contrast to the previous models, Model E does not include 
a contribution from the Comptonization of soft photons (i.e. 
the COMPTT model). Alternatively, the excess seen at soft 
energies is modelled using a blurred reflector in addition 
to reflection consistent with distant emission. The ionized 
reflection model reflionx is convolved with a kerrdisk 
kernel, modelling a relativistically blurred reflection spec- 
trum from the accretion disc originating from very near to 
the black hole instead of the kerrdisk model i.e. the ker- 
RCONV model (Brenneman & Reynolds 2006) . The reflection 
components consist of: (reflionx + kerrconv*reflionx). 
The assumptions in using this convolution model are con- 
sistent with those used in Model D, the inner radius of the 
accretion is assumed to extend down to risco and outer ra- 
dius fixed at 400 r ms for a disc of uniform emissivity. Narrow 
emission lines from ionized matter are also included where 
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Table 7. Fit parameters from Model D to Suzaku XIS, HXD and BAT data from Swift. Line energies are quoted in the rest frame. 
* denotes a frozen parameter. a POWERLAW normalization given in units (10 -2 phkeV -1 cm -2 s — x ). b REFLIONX normalisation given in 
units 10 -5 . 
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required, while the narrow 6.4 keV core is accounted for by 
the second unblurred reflionx component. 

Model E is a significantly worse fit in all 4 AGN with a 
soft excess (by an average of A% 2 ~ 23 compared to Model 
D), requiring a significant amount of blurring to account for 
the excess at low energies (see Table [S]). Given the amount 
of blurring required, parameters such as the emissivity index 
and the spin parameter approach more extreme values for all 
objects with a soft excess (e.g. q > 4 and a > 0.9). Whilst 
Model E is a good fit to the data for Mrk 335 and NGC 
7469, Model D provides a better fit and without the use of 
such extreme parameters. 

In addition, the best-fit obtained with Model E for Ark 
120 requires an inclination of i — 48°, whilst not unreason- 
able, differs from the inclinations measured using Models C 
& D. However freezing the inclination at i = 35° worsens the 
fit further ( X l = 824.6/648 compared to xl = 770.2/647). A 
similar dual-reflector fit as used here is considered by Nar- 
dini et al. (2010) in which the emissivity index is frozen at 
q — 5.0 compared to q = 4.1+g'g measured here. Nardini et 
al. (2010) also note a relatively high inclination of i° — 57+^ 2 
obtained in their analysis is likely too large for an object 
such as Ark 120 although it possible for the accretion disc 
and obscuring material to be misaligned. One area of dis- 



agreement between the results obtained with Model E and 
the Nardini et al. (2010) fit is the determination of the spin 
parameter. Nardini et al. finding 0.24 < a < 0.93 whereas 
here only a lower limit of a > 0.97 is found, this may be due 
to the degeneracies between the spin and emissivity index 
which are both allowed to vary in Model E. 

Similarly, NGC 7469 is best fit with a kerrconv incli- 
nation parameter of i° — 70^3 whereas Models C & D sug- 
gest that the disc is inclined at i ~ 23° to the observer, wors- 
ening the fit by Ax 2 ~ 80. These high inclinations are likely 
to be driven by the need to model the relatively smooth soft 
excess. 

For the case of MCG-02- 14-009, Model E provides an 
approximately equal quality of fit compared to Models C & 
D. No significant amount of blurring is required to model 
the spectrum, obtaining parameters similar to those pre- 
viously i.e. low emissivity and an unconstrained spin pa- 
rameter. The narrow FeXXVI ionized emission line is also 
present, in line with findings from Models C & D. The 
SWIFT J2127.4+5654 spectrum is fitted very well with 
Model E (xt = 830.4/865). Extreme parameters are not re- 
quired to fit the data, which is in agreement with results 
from the previous disc line models. This may be due to 
the lack of any notable soft excess in these two objects. 
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Figure 5. uFv plots of Model D indicating the strength of any soft excess and relativistic line emission from the kerrdisk models. 
Consisting of POWERLAW + COMPTT (where required) + unblurred REFLIONX + KERRDISK + narrow ionized lines as required. XIS data 
is in black, HXD in red and BAT data is represented by blue circles. 



Since a significantly blurred spectrum simulates an excess 
at lower energies, a high emissivity and near maximally ro- 
tating black hole would provide an over-excess at soft ener- 
gies, inappropriate for objects such as MCG-02- 14-009 and 
SWIFT J2127.4+5654. Nonetheless Miniutti et al. (2009) 
find q = 5.3±i;I in SWIFT J2127.4+5654 compared to 
q — 2.2^0 4 here, this may be due to the high level of Galac- 
tic absorption and the additional small amount of intrinsic 
absorption at the redshift of the source. An increase in the 
emissivity index could be compensated for by increased ab- 
sorption which is found in the Miniutti et al. (2009) analysis. 



3.6 Model F - Blurred reflector (2) 

Model F starts with Model B as the base model but instead 
blurring the single reflionx reflection spectrum with the 
kerrconv convolution model to model any broad residu- 
als present in the spectra. The FeK line complex can still 
be modelled in this way since the reflionx model includes 
Fe Ka emission and blurring the spectrum emulates the re- 
sulting profile from a kerrdisk model. Since the narrow 
Fe Ka core included within reflionx is now relativistically 
blurred, a narrow Gaussian of fixed width 10 eV was added 
to ensure that the narrow 6.4 keV is still modelled. This may 
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Table 8. Fit parameters from Model E to Suzaku XIS, HXD and BAT data from Swift. (A) represents the unblurred REFLIONX and 
(B) represents the blurred REFLIONX. The ionization parameter £ is given in units erg cms -1 . a POWERLAW normalization given in units 
(10 — 2 phkeV -1 cm -2 s _1 ). b REFLIONX normalization in units 10~ 5 . * denotes a parameter frozen at the best-fitting value from Model 
D. In some cases the spin parameter a could not be constrained, denoted by - . 
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represent the case where the 6.4 keV line is observed from 
Compton-thin matter, such as the BLR or NLR. Contrary to 
Model E, the soft excess (where present) is modelled using 
the COMPTT soft photon Comptonization model. 

Any ionized emission due to FeXXV and FeXXVI re- 
quired in Model D is also found to be required in this model. 
In general a good fit to all objects is obtained, with Model F 
clearly providing a better fit to the data compared to Model 
E for Fairall 9, NGC 7469 and SWIFT J2127.4+5654. The 
kerrconv parameters are also consistent within error bars 
with the emissivity index, inclination and spin parameter 
found using the kerrdisk line profile previously, yielding 
typically slightly lower emissivity indicies. 

Note that the best-fitting spin parameter value with 
Model F for Fairall 9 is a = 0.401q |p (quoted at the 75% 
confidence level and is unconstrained at the 90% confidence 
level) in agreement with a = 0.44l8' l j'i found in Model D 
(at the 90% confidence level) and with a = O.5I8 3 found by 
Schmoll et al. (2009). 

Given the independent modelling of the soft excess 
within this model, it is interesting to note that a = 0.72l8:i 7 
is obtained for NGC 7469 at the 90% confidence level. This 
value of the spin parameter is in agreement with that found 
in Model D, although providing a slightly worse constraint. 



3.7 Comparison with XMM-Newton results 

Analysis of the data for these objects obtained with XMM- 
Newton (all objects other than SWIFT J2127.4+5654, see 
Table [2]) yields no additional constraints upon the properties 
of the Fe K region, the reflection component is also difficult 
to constrain due to the lack of hard X-ray data. As with 
the Suzaku data, the FeXXVI narrow ionized emission line 
is found throughout most objects and the FeXXV line is 
relatively uncommon, observed only in Ark 120 and Fairall 
9. However, the results obtained with a kerrdisk relativistic 
line profile are consistent with those obtained with Model 
D in the Suzaku data. Similarly, a good fit is obtained with 
this model for all of the objects in the sample, with maximal 
spin ruled out from the fits. The 0.5-10.0 keV spectrum is 
generally well described by a combination of narrow distant 
emission and a small broad relativistic component, typically 
EW < 150 eV. 

In particular, the XMM data for MCG-02-14-009 yields 
very similar results to the Suzaku data. The main differ- 
ence being that a statistically significant FeXXVI line is 
not present in the XMM-Newton data. All measurements 
obtained with a kerrdisk profile are consistent at the 90% 
confidence level with the Suzaku data. Both data sets there- 
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Table 9. Fit parameters from Model D to XMM-Newton EPIC-pn data. * denotes a frozen parameter. In some cases the spin parameter 
a could not be constrained, denoted by — . a POWERLAW normalization given in units (10 — 2 phkeV - 1 cm -2 s" 1 ). b Flux for COMPTT 
quoted over the 0.5-10.0keV range in units 10 —11 erg cm -2 s — 1 . c REFLIONX normalisation given in units 10 — 5 . 
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fore reveal a relatively simple 0.5-10.0 keV spectrum, there 
is no observed soft excess, it is well described by distant re- 
flection in addition to a moderately broad relativistic line 
profile from neutral Fe occuring at r- m > 2.5 r g (obtained 
from the application of the LAOR model) with emissivity 
q = 2.8^0.8 > as seen in Tabled 

These same observations of Ark 120 and Fairall 9 have 
previously been modelled with the blurred reflector ap- 
proach, consisting of a combination of narrow distant emis- 
sion and a smeared ionized reflector under the assumption 
of a Schwarzschild geometry, see Brenneman & Reynolds 
(2009). Consistent results are found for Fairall 9 for which 
narrow ionized emission lines are also found, however the 
authors were unable to constrain the emissivity index, mea- 
sured at q < 1.8 here and the spin value requires that emis- 
sion occurs at r- m > 1.8 r g . The results found here for Ark 
120 suggest consistent measurements for the emissivity in- 
dex with Brenneman & Reynolds (2009). In comparison with 
the Suzaiu data for this object, both approaches are consis- 
tent with Model D. Ionized emission lines are also found in 
Ark 120 throughout both analyses. The inner radius of emis- 
sion is generally consistent, occurring at n n < 80 r g from 
Brenneman & Reynolds (2009) and r- m > 2.5 r g as mea- 



sured here from the spin parameter within Model D (Table 

[5}. 



4 DISCUSSION 

This small sample of six AGN includes typically bare Seyfert 
1 galaxies featuring little or no intrinsic absorption. This 
property of these AGN is important since the spectra of 
these objects is simpler to model without complicating fac- 
tors (such as absorption) allowing the observer to draw con- 
clusions about the fundamentals of accretion disc properties 
and basic features of the Fe K region. These conclusions will 
therefore be less model dependent up on how the warm ab- 
sorber is modelled (Turner & Miller 2009). It is important 
to effectively model the spectra of these AGN and assess 
the likely origin of various components of the spectrum be- 
fore we can proceed to draw conclusions about more com- 
plicated AGN. Modelling of the broad-band continuum is 
also essential prior to analysis of the Fe K region, this has 
been achieved using data from both HXD and BAT hard 
X-ray detectors onboard Suzaku and Swift respectively. Al- 
lowing the spectrum spanning 0.5-100.0 keV to be modelled 
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Figure 6. i/Fu plots of Models E & F for Mrk 335 showing to degree to which relativistic blurring of the reflection component is required 
according to differing interpretations of the soft excess. Note that the significant blurring required in Model E to account for the soft 
excess reduces the accuracy of the fit to the Fe K region features. XIS data is in black, HXD in red and BAT is represented by blue 
circles. 



with better constraints upon distant reflection components 
modelled by pexrav and reflionx, whilst maintaining ef- 
fective modelling of the soft excess. This approach furthers 
that taken by Nandra et al. (2007) in which EPIC-pn data 
from XMM-Newton is analysed in the 2.5-10.0 keV range. 

The presence of narrow ionized emission lines due to 
Fe XXV and Fe XXVI assessed prior to modelling any broad 
residuals in the Fe K region has an important effect upon the 
parameters obtained with broad disc line profiles. In some 
spectra neglecting to model these lines (if present) can ac- 
centuate any apparent broad Fe K line, particularly those re- 
sulting from ionized species of Fe. Most commonly occurring 
in the objects analysed here is the 6.97 keV line from H-like 
Fe which is observed in all objects, except NGC 7469 and 
SWIFT J2127.4+5654. This is a surprising result given the 
rarity of such lines found by Nandra et al. (2007) in which 
only 2/26 sources showed evidence for these lines (although 
more common in work by Bianchi et al. 2009). This may be 
due to higher quality (longer exposure) data obtained with 
Suzaku allowing these lines to be more easily distinguished 
from broad residuals at these energies. Similarly, neglecting 
to include a narrow 6.97 keV emission line where present 
forces the relativistic line profile to have the emissivity in- 
dex increased to particularly large values (q > 6) and the 
inclination is slightly increased to ~ 45° — 55° whereby the 
blue wing of the profile is forced to model the narrow excess 
at ~6.97keV. 

In agreement with Nandra et al. (2007), however, is the 
rarity of the 6.7 keV line from He- like Fe which only fea- 
tures significantly here in Mrk 335. Models A & B earlier 
suggested the possible presence of FeXXV in Fairall 9 and 
SWIFT J2127.4+5654, however it was found that excesses 
at these energies coincided with the peak in the blue-wing 
of the relativistic line profile. This suggests that high (i.e. 
calorimeter) resolution data is necessary to determine the 
presence of this particular emission line particularly when 
employing a relativistic line profile, such as that which would 
be obtained with Astro-H (see Kelley et al. 2010). 

The introduction of the kerrdisk model over a model 
considering emission purely from distant material (i.e. 
Model B) improves the fit to the six objects by an average 



of Ax 2 ~ 20. This implies that broad residuals are a sta- 
tistically significant feature in all of these objects (although 
less so in MCG-02- 14-009). 

4.1 Typical parameters of the accretion disc 

The average inclination of the accretion disc to the observer 
as inferred by the kerrdisk line profiles in Model D is i — 
33° ± 4°, this is consistent with Nandra et al. (2007) who 
find i — 38° ± 6°. Also from Model D we find an average 
emissivity index of q = 2.3 ± 0.2, much lower than those 
used in previous work for some of these objects e.g. Fairall 9 
(Schmoll et al. 2009) and SWIFT J2127.4+5654 (Miniutti et 
al. 2009), particularly when using blurred reflection models 
to model the whole continuum. The line profile produced 
within Model D also suggests a broad line profile with an 
average equivalent width EW = 119± 19 eV consistent with 
an average EW = 91.3 ± 12.8 found by Nandra et al. (2007) 
using a Gaussian to model the broad residuals. Whilst here 
for Models C & D the centroid rest energy of the line profile 
is not allowed to extend below 6.4 keV in the rest frame, the 
basic parametrization within Model A suggests an average 
broad component line energy of 6.29 ±0.03 keV compared to 
E Ka = 6.27 ± 0.07keV according to Nandra et al. (2007). 

The low average emissivity index of the objects in this 
small sample may be due to a number of factors: indepen- 
dent modelling of the soft excess through a Comptonization 
of soft photons; modelling of ionized emission lines where 
present in the data and the assumption that emission ex- 
tends down to nsco (the kerrdisk model). In accordance 
with this, within Models E & F facilitating the use of a 
blurred reflector, the average emissivity index is also de- 
pendent upon the way in which the soft excess is mod- 
elled. Model F (in which the soft excess is modelled with 
a COMPTT component, as in Models A-D) suggests an av- 
erage q = 2.0 ±0.1 in general agreement with Model D. 
However Model E (in which purely the blurred reflection 
component is tasked with modelling the soft excess) sug- 
gests an average q — 5.0 ± 0.7 for those objects featuring a 
soft excess and an average q = 2.0 ± 0.4 for those without, 
namely MCG-02-14-009 and SWIFT J2127.4±5654. Thus 
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Figure 7. Confidence plots from Model D for the spin parameter a for those objects in which it could be constrained. Dashed lines 
representing 90% (black), 95% (red), 99% (green) and 99.5% (blue) confidence levels. Note that the emissivity index q is a free parameter. 



the high emissivity in Model E for some objects is driven by 
the need to fit a featureless soft excess e.g. see Figure [6] A 
blurred reflection model with high emissivity therefore ap- 
pears to be ruled out given Model E produces a significantly 
worse fit in all 4 AGN with a soft excess. 

Model C alternatively models reflection off a relativistic 
accretion disc under the assumption of a maximally rotating 
central black hole, without the inner radius of emission fixed 
at the ISCO. This approach also suggests a similarly rela- 
tively small equivalent width of EW = 102 ± 14 eV and an 
average emissivity q — 3.2 ± 0.4 (for those objects for which 
the emissivity could be constrained), higher than that de- 
rived from Model D and likely due to the fact that is 
not fixed at nsco- The average inclination of the accretion 
disc is also very similar to Model D above with i — 34° ± 3° . 
Given the assumptions made within this model, the inner 
radius of emission is found to originate at tens of r s rather 
than < 6r g at an average n n = 39±8 r g . Indeed, in Ark 120, 
Fairall 9 and MCG-02- 14-009 Model C and this interpreta- 
tion of the accretion disc and SMBH provides a marginally 
better fit. However given the degeneracies between q, a and 
ri n it would be impossible to reasonably constrain any of 
these values without assumptions similar to those mentioned 
above. 



4.2 Model dependent Black Hole spin 

The results for determining the spin of the central black 
holes of the AGN in this sample suggest that the spin de- 
rived is very much dependent upon which interpretation of 
the FeK line region is followed. Modelling the soft excess 



through an independent model such as compTT tends to 
yield low to intermediate spin constraints for the objects 
in this sample, the exceptions being Ark 120 and MCG-02- 
14-009 in which only an upper limit could be placed. The 
employment of this interpretation also yields low to interme- 
diate emissivity indicies for the accretion discs (i.e. q ~ 2). 
Modelling the 0.5-100.0 keV spectrum with a blurred reflec- 
tion component and no other modelling of the soft excess 
(i.e. Model E) gives particularly high values of both emis- 
sivity index and spin parameter, typically q > 4 and a > 0.9, 
but these fits are statistically ruled out here. 

This is, however, only in objects featuring a soft ex- 
cess. In MCG-02-14-009, which has no obvious excess over 
a powerlaw at low energies, the derived parameters are very 
similar throughout Models D, E & F. Similarly for SWIFT 
J2127. 4+5654, the parameters obtained with these three 
models are all consistent and indeed the resulting spin pa- 
rameter is also consistent with previous findings by Miniutti 
et al. (2009) although they find a much higher emissivity in- 
dex is required (q — 5.3^\'\ compared to q = 2.2j^'g from 
Model E). 

As discussed previously, the spin constraint obtained 
here for Fairall 9 agrees with that found by Schmoll et al. 
(2009) , but only in the case where they ignore the spectrum 
below 2keV (a = 0.5±g.l compared to a = 0.441^* found 
here in Model D). According to Model D, an intermediate 
spin of the central black hole within NGC 7469 is found, 
a = 0.691q'q9- This is also consistent with that found within 
Model Fa — 0.72+q'^, suggesting that the spin of this ob- 
ject is indeed a ~ 0.7. 

These findings therefore suggest that the two differing 
interpretations give concurring spin constraints only when 
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the blurred reflector is not required to model the soft excess, 
leading to the conclusion that the component responsible for 
the soft excess must be independently fitted before conclu- 
sive spin constraints can be made. Even so, the origin of 
the soft excess is not currently known, e.g. see Gierlinski & 
Done (2004). The relatively constant temperature of the soft 
excess versus BH mass suggests that it may not arise from 
direct thermal emission from the accretion disc, since the ac- 
cretion temperature properties should scale with M BI J' 4 in a 
standard accretion disc. An atomic origin of the soft excess 
has been suggested, however no obvious spectral features are 
seen in high resolution data meaning that if atomic emission 
is responsible it must be significantly relativistically blurred 
such as here in Model E (Ross, Fabian & Ballantyne 2002; 
Fiore, Matt & Nicastro 1997). Alternatively, it is suggested 
that the soft excess could arise from relativistically blurred 
absorption (Gierlinski & Done 2004) from a differentially ro- 
tating and outflowing disc wind (Murray & Chiang 1997). 
The smeared absorption creates a smooth hole in the spec- 
trum resulting in the apparent soft excess and a hardening 
at high energies. Whilst in this paper we do not aim to 
determine the origin of the soft excess, we find that a sim- 
ple parametrization of the soft excess continuum through a 
model such as compTT provides a better fit to the spec- 
tra than an atomic origin from a highly blurred reflection 
component. 

4.3 Reflection off a disk wind? 

One possible alternate origin of the observed broad compo- 
nent in AGN spectra is through reflection of primary con- 
tinuum photons off a Compton thick disc wind, producing 
emission, absorption and broad features in the FeK band 
(Sim et al. 2010). Fast outflows can add significant com- 
plexity to the Fe K region and can reproduce a wide range 
of spectral signatures owing to the many differing physical 
conditions. Sim et al. (2010) note that such outflows can, 
however, have a significant affect upon the soft X-ray spec- 
trum resulting in features such as highly blueshifted absorp- 
tion lines, although these absorption are weaker than the 
features seen in the FeK band. The lines from a disc wind 
may not be observed in objects such as those in this sample 
if they are observed relatively face-on, however even if we are 
not looking directly through the line of sight to the wind, it 
may be possible to observe broadened emission features in 
reflection. For instance Sim et al. (2010) predict broadened 
FeK emission down to ~ 5keV in their disc wind spectra, 
with equivalent widths typically several tens of eV. 

5 CONCLUSIONS 

Resulting from the work above on this sample of six 'bare' 
Seyfert AGN, we conclude the following: 

(i) Narrow emission from distant material is very impor- 
tant when modelling the FeK region and neglecting to in- 
clude these components where present can have a significant 
affect upon the accretion disc parameters obtained with sub- 
sequent models. The narrow neutral 6.4 keV core is ubiqui- 
tous amongst these six AGN. 

(ii) Ionized emission at ~ 6.97 keV due to FeXXVI is 



more common in these spectra than previously thought (e.g. 
Bianchi et al. 2004 and Nandra et al. 2007) although Bianchi 
ct al. (2009) find them more common. It is present here in 
4/6 objects in the Suzaiu data and it is also present in the 
same objects during the 5 XMM-Newton observations in- 
cluded here. 

(iii) Emission at ~ 6.70 keV due to FeXXV is still rela- 
tively uncommon amongst these AGN, present in only 1/6 
SuzaJcu and 2/5 XMM-Newton observations. However, it is 
likely that this emission line may be more common amongst 
AGN since this is also the energy at which the blue-wing 
from a relativistically blurred neutral 6.4 keV line occurs 
with typical model parameters (e.g. i = 30° —40°, q = 2 — 3), 
making it difficult in most cases to distinguish between these 
two components in CCD resolution spectra. 

(iv) The fit to all objects in the sample is significantly 
improved with the introduction of a relativistic line profile, 
such as that from a kerrdisk model. However this broad 
emission is weaker after effective modelling of narrow com- 
ponents (mean EW = 119 ± 19 eV). 

(v) The average emissivity index of the accretion disc is 
a low to moderate q = 2.3 ± 0.2 when using compTT to 
model the soft excess. The emissivity index typically scales 
as R~ 2 3 , consistent with an RT 2 law for the illuminating 
continuum and therefore implying that strong GR effects 
(such as light bending) may not be required in these ob- 
jects. The assumption that emission extends down to nsco 
is still valid, however, since the flat emissivity of the disc 
indicates that emission is not centrally concentrated. This 
is consistent with emission occuring a typically tens of r g 
rather than within < 6 g . 

(vi) It is essential to effectively model the continuum and 
the soft excess prior to attempting to constrain parameters 
of the accretion disc and inner regions of AGN. Independent 
modelling of the soft excess with models such as compTT 
or a second soft powerlaw, rather than blurring reflection 
from inner parts of the accretion disc as a way to model the 
soft excess yields consistent results with respect to AGN 
with and without an excess at lower energies. 

(vii) Purely using a blurred reflection component to 
model the soft excess results in typically higher values of 
spin and emissivity index. This would lead to a skewed dis- 
tribution amongst AGN towards the higher values of these 
parameters for those objects with a soft excess i.e. a > 0.9 
and q > 4. This would indicate that the regions of the accre- 
tion disc responsible for the soft excess (in an atomic origin 
of the soft excess) are not the same as those responsible 
for features in the FeK region. In this small sample of six 
AGN, the fit is significantly worse using this interpretation 
for those objects with a soft excess, namely Ark 120, Fairall 
9, Mrk 335 and NGC 7469. 

(viii) New constraints suggest further intermediate spin 
values for Mrk 335 and NGC 7469 of a = 0.70to.ln and 
a = 0.69±g;g| respectively. The spin value for SWIFT 
J2127. 4+5654 found here is consistent with that found by 
Miniutti et al. (2009) , whereas the spin found here for Fairall 
9 is consistent with the Schmoll et al. (2009) analysis above 
2keV. Only upper bounds for the spin parameter can be 
placed for Ark 120 and MCG-02-14-009, requiring that emis- 
sion does not occur within 2.02 r g and 2.45 r g respectively. 

(ix) Overall zero spin cannot be ruled out at a particu- 
larly high confidence level in all objects, as can be seen in 
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Table [7] None of the objects in the sample show overwhelm- 
ing evidence for a Kerr geometry and significantly centrally 
concentrated emission, maximal spin is ruled out in most 
cases (Models D & F). 
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